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METHODS  FOR  COMPUTING  THE  PATTERN 
OF  COMPRESSIBLE  FLUID  FLOW  IN 
NOZZLES  OR  PAST  BODIES 

FINAL  REPORT 

ON 

AFOSR  CONTRACT  NO.  AF-F44 62 0 -74 - C- 0062 

Cathleen  S.  Morawetz  , 

A  paper  on  using  the  method  of  relaxation  for  a  tran¬ 
sonic  type  mixed  equation  corresponding  to  transitions  from  sub¬ 
sonic  to  supersonic  flows- yas  completed  [5].  It  was  shown  how 

vr  ‘ 

to  apply  this  technique^  to  obtain  convergent  iteration  schemes 
for  an  equation  of  Tricomi  type  in  a  rectangle. 

Amiram  Harten 

'Research  on  the  method  of  artificial  compression  con¬ 
sisted  of  two  parts:  . 

(1)  Development  and  analysis  of  the  scheme; 

(2)  Applications  to  computation  fluid  dynamics  and, 
in  particular,  to  flows  with  chemically  reacting 
fluids . 

(1)  The  development  and  analysis  part  of  the  study  is 
summarized  in  a  series  of  four  articles.  The  first  installment 
of  this  series  [6]  deals  with  artificial  compression  for  single 
conservation  laws,  and  has  appeared  in  the  Communications  on  Pure 
and  Applied  Mathematics.  In  the  second  installment  of  this  series 
he  discusses  artificial  compression  for  one-dimensional  systems  of 
conservation  laws.  This  method  is  motivated  by  Liu's  [1]  generali¬ 
zation  of  Oleinik's  entropy  condition.  This  article  appeared  in 
the  Communications  on  Pure  and  Applied  Mathematics.  The  third 
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installment  of  the  series  [6]  describes  a  technique  to  apply  arti¬ 
ficial  compression  in  conjunction  with  high  order  accurate  schemes. 

A  study  of  local  monotonicity  of  finite  difference  schemes  at  dis¬ 
continuities  yields  a  necessary  condition.  This  condition  is  used 
to  construct  a  hybrid  scheme  in  which  a  high  order  scheme  is  auto¬ 
matically  switched  into  a  monotone  scheme  at  discontinuities. 
Numerical  results  presented  in  this  paper,  as  well  as  those  in  a 
comparative  study  by  Gary  Sod  [2],  demonstrate  a  remarkable  improve¬ 
ment  in  the  performance  of  the  corrected  schemes.  This  paper  has 
appeared  in  the  Mathematics  of  Computation  [7].  In.  the  fourth  and 
last  installment  of  this  series,  artificial  compression  is  analyzed 
in  the  multi -dimens ional  case,  and  its  performance  demonstrated  on 
realistic  problems  such  as  supersonic  flows  around  cylindrical  bell 
shaped  bodies  and  interaction  of  an  impinging  shock  with  a  bow 
shock. 


S.  H.  Burstein  and  A.  Harten 


Burstein  and  Harten  made  a  study  of  the  feasibility  of 
replacing  shock  fitting  techniques  for  computations  of  chemically 
reacting  fluids  by  (i)  artificial  compression  type  methods  and 
(ii)  image  enhancement  methods.  When  computing  chemically  reacting 
flows  it  is  extremely  important  to  have  a  monotone  and  sharp  front. 
The  work  in  this  direction  consisted  of  two  main  approaches: 

(a)  An  attempt  to  apply  artificial  compression  directly  as  before. 
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compTession  method  by  substituting  stTict  conservation  by  conser¬ 
vation  in  the  mean.  This  work  is  motivated  by  the  random  choice 
method  of  Chorin  [4];  (c)  Transforming  the  differential  equations 
using  Tschebyschef f  polynomials  of  the  first  kind  T  (x)  ,  as  a 
representation  for  the  flux  terms.  The  mesh  spacing  is  chosen  to 
follow  a  cosinusoidal  variation  so  that  T  (5)  cos  £  which 
allows  for  the  use  of  the  Fast  Fourier  Transform.  The  advantage 
of  the  method  when  applied  to  cylindrically  expanding  shocked 
flows,  is  twofold:  first,  mesh  spacing  is  dense  behind  the  shock 
where  chemical  reaction  takes  place  and  at  the  boundary  (which 
drives  the  rarefaction)  and  second,  since  the  computation  of  the 
divergence  of  the  flux  takes  place  in  wave  number  space,  special 
filtering  methods  used  in  image  enhancement  processes  can  be 
tested  for  applicability  to  the  moving  shock  problem. 

Contacts  with  Kirkland  Base  were  maintained  and  all 
material  for  incorporating  the  artificial  compression  method  into 
existing  codes  was  sent  there. 
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Heinz-Otto  Kreiss 

B.  Gustafsson  and  Heinz-Otto  Kreiss  have  written 
the  paper.  Boundary  conditions  for  time-dependent  problems 
with  an  artificial  boundary.  This  is  joint  work  with  B. 
Gustafsson  to  appear  in  the  Journal  of  Computational  Physics. 
Some  of  the  more  interesting  results  are: 

1.  One  should  not  overspecify  boundary  conditions. 

2.  Upstream  differencing  can  give  completely 
wrong  results.  The  numerical  solution  can 
appear  smooth  but  have  nothing  to  do  with 
reality . 

3.  To  avoid  these  troubles  the  authors  have 
derived  asymptotic  expansions  which  are 
connected  with  the  difference  schemes. 


Consider  the  potential  equation  for  potential  flow: 


(1  '  *  a  -  pr>V 


+  -  =  0 


and  write  it  as  a  first  order  system.  To  improve  the  efficiency 
of  computational  nozzle  flow  it  would  be  very  desirable  to 
approximate  by  the  "box  scheme,"  which  is  the  most  compact  second 
order  accurate  approximation.  Unfortunately,  numerical  results 
show  that  the  straight  forward  application  of  the  method  does 
not  work.  Work  has  started  to  explain  the  difficulty  and  to 
modify  the  method. 


Michael  Mock 

/ 

Michael  Mock  studied  singular  non-linear  elliptic 
problems  by  the  finite  element  method.  In  studying  the 
strength  of  elastic  structures  a  great  deal  of  use  has  been 
made  of  the  finite  element  method  to  approximate  the  partial 
differential  equations  which  govern  the  equilibria.  Mock 
attacked  the  problem  with  a  view  to  making  use  of  his  long 
experience  with  semi-conductors  which  have  several  common 
mathematical  traits  with  elastic  structures,  notably  very 
sharp  regions  of  transition,  [8],  [9]. 
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